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Abstract 

Stratospheric intrusions have been the interest of decades of research for their ability to 
bring stratospheric ozone (O3) into the troposphere with the potential to enhance sur- 
face Og concentrations. However, these intrusions have been misrepresented in models 
and reanalyses until recently, as the features of a stratospheric intrusion are best iden- 
tified in horizontal resolutions of 50 km or smaller. NASA’s Modern-Era Retrospective 
Analysis for Research and Applications Version-2 (MERRA-2) reanalysis is a publicly- 
available high-resolution dataset (~50 km) with assimilated O3 that characterizes O3 
on the same spatiotemporal resolution as the meteorology. We demonstrate the science 
capabilities of the MERRA-2 reanalysis when applied to the evaluation of stratospheric 
intrusions that impact surface air quality. This is demonstrated through a case study 
analysis of stratospheric intrusion-influenced O3 exceedences in spring 2012 in Colorado, 
using a combination of observations, the MERRA-2 reanalysis and the Goddard Earth 
Observing System Model, Version 5 (GEOS-5) simulations. 


1 Introduction 


Surface ozone (O3) is harmful to human health and agriculture [Scherrer et al., 2006; 
Krzyzanowski and Cohen, 2008]. Near the surface, O3 is termed a secondary pollutant 
since it is a product of the photochemical reaction with precursors such as nitrogen ox- 
ides (NO,; NO and NO), carbon monoxide (CO) and non-methane hydrocarbons which 
have both man-made and natural emission sources in the troposphere. Therefore, in or- 
der to reduce near surface Og concentrations, communities must reduce anthropogenic 
pollution sources. However, the injection of stratospheric O3 into the troposphere, known 
as a stratospheric intrusion (SI), can also lead to concentrations of ground-level O3 ex- 
ceeding the national ambient air quality standard (NAAQS) set by the Environmental 
Protection Agency (EPA), especially at high elevations [e.g., Langford et al., 2009, 2015; 
Lin et al., 2012, 2014; Yates et al., 2013; Zhang et al., 2014]. In October 2015, the EPA 
revised the US NAAQS for daily maximum 8-hour average (MDA8) O3 from 75 parts 
per billion by volume (ppbv) to 70 ppbv [U.S. Environmental Protection Agency, 2015]. 
Therefore, it is crucial that we are able to understand, model, and predict SIs and their 
potential impact on surface O3 concentrations. 


SIs form as a result of the tropopause being drawn down below the jet stream, re- 
ferred to as tropopause folding, often associated with an upper-level trough. SIs are char- 
acterized by O3-rich [e.g., Danielsen, 1968; Shapiro, 1974, 1980; Holton et al., 1995; Brown- 
ing, 1997] and CO-poor [Fischer et al., 2000] air, with relatively high levels of potential 
vorticity (PV) [Holton et al., 1995] and low levels of water vapor often observed in satel- 
lite imagery as a “dry slot” [e.g., Bader et al., 1995; Wimmers et al., 2003]. Therefore, 
tropopause folds can lead to the mixing of stratospheric and tropospheric air with dif- 
ferent chemical and meteorological properties at low altitudes [e.g., Danielsen, 1980; Shapiro, 
1980; Holton et al., 1995], remaining behind a mid-latitude cyclone’s surface cold front 
Browning, 1997; Bethan et al., 1998; Cooper et al., 2001; Knowland et al., 2015]. The 
west coast of the USA is located at the end of the Pacific Ocean storm track [e.g. Hoskins 
and Hodges, 2002], a region favorable for stratosphere-to-troposphere transport of O3 
James et al., 2003; Sprenger and Wernli, 2003; Stohl et al., 2003; Skerlak et al., 2014]. 

On the lee-side of the Rocky Mountains, cyclones form (or redevelop) supported by upper- 
level troughs [McClain, 1960; Carlson, 1991]. However, the descending motion associ- 

ated with the upper-level trough can still be a strong feature in the troposphere over the 
Rocky Mountains, prior to the identification of a surface low pressure system. In the upper- 
level flow, the troughs can form closed lows and even become “cut-off” from the west- 

erly flow [Palmén and Newton, 1969]. This can result in the prolonged influence of the 
tropopause folds on tropospheric O3 concentrations over a region [Lin et al., 2012; Yates 

et al., 2013] until the cut-off low (COL) dissipates or is reabsorbed into the mean flow 
[Nieto et al., 2008]. During the winter and spring, there is a build-up of O3 in the lower 


stratosphere, and this leads to SIs having the largest influence on surface O3 in the spring 
[Danielsen and Mohnen, 1977; Holton et al., 1995; Monks, 2000]. 


For over 40 years, studies have observed the injection of O3-rich air into the tro- 
posphere within tropopause folds over western USA [e.g., Lovill, 1970; Shapiro, 1980; 
Langford et al., 1996, 2009, 2012; Wimmers et al., 2003; Cooper et al., 2004; Lefohn et al., 
2011] with recent studies focusing on the impact of SIs on O3 air quality exceedences in 
the high elevation communities of the Rocky Mountains [e.g., Langford et al., 2009, 2015; 
Lin et al., 2012, 2014; Yates et al., 2013; Zhang et al., 2014]. Langford et al. [2009] fo- 
cused on the transport of stratospheric O3 in Colorado’s Front Range during the spring 
of 1999 using lidar and surface measurements of O3. They identified high concentrations 
of O3 in the mid-troposphere down to the surface as a result of tropopause folds asso- 
ciated with upper-level troughs in the region. Lin et al. [2012] utilized the abundance 
of vertical observations from ozonesondes and lidar taken during the 2010 NOAA Cal- 
Nex field campaign in California as well as ground-based measurements throughout west- 
ern USA, in conjunction with a model study to quantify the stratospheric fraction of air 
that impacts NAAQS exceedence events. Using the NOAA Geophysical Research Lab- 
oratory (GFDL) Atmosphere Model version 3 (AM3) with fully coupled stratosphere- 
troposphere chemistry at ~50 km resolution, Lin et al. [2012] attributed 50-60 % of to- 
tal modelled surface O3 in spring 2010 (as much as 20-40 ppbv of additional O3 during 
4 deep intrusions) to stratospheric origins on exceedence days. Using a coarser resolu- 
tion model (~200 km), Lin et al. [2015] extended the analysis to April and May during 
a 23-year period (1990-2012) and found the average stratospheric O3 contribution is 15— 
25 ppbv of western US surface O3. 


While the impact of SIs on surface O3 in the western US is well documented, sim- 
ulating and predicting such events remains challenging. The resolutions of current global 
meteorological analyses (~10-50 km) are sufficient for resolving the dynamical evolution 
of SIs, however these models typically contain very limited representations of trace gases 
like Og. Reanalyses have been used in numerous studies to explore the frequency, spa- 
tial variations and structure of SIs [e.g., Stohl and Trickl, 1999; Waugh and Polvani, 2000; 
Sprenger and Wernli, 2003; Lefohn et al., 2011; Reutter et al., 2015; Nath et al., 2016], 
however, there are very few such studies which also use reanalysis O3 [Skerlak et al., 2014; 
Zanis et al., 2014; Knowland et al., 2015, 2017; Ott et al., 2016; Ryoo et al., 2017]. It 
is our objective to investigate whether NASA’s Modern-Era Retrospective Analysis for 
Research and Applications Version-2 (MERRA-2) reanalysis, which is similar to NASA’s 
Global Modeling and Assimilation Office (GMAO) operational forecasting system, is able 
to capture the dynamical features of a SI, in particular the isentropic descent of elevated 
Os within and below the tropopause fold. Such datasets would support air quality agen- 
cies for more rapid identification of the impact of stratospheric air on ground-level O3 
[Kaldunski et al., 2017] separate from local sources or the long-range transport of O3 [Ryoo 
et al., 2017]. The focus of this study will be on springtime (March - June (MAMJ)) O3 
air quality exceedences in 2012 which were identified by the EPA as having direct con- 
nection with SIs [US EPA AQS database, 2017]. 


2 Data 
2.1 Observational datasets 


In the spring of 2012, there were seven days when the MDA8 O3 [EPA AirData, 
2016] at the Rocky Mountain National Park (RMNP) Long’s Peak monitoring station 
(40.27°N, 105.54°W, 2742 m, Air Quality System (AQS) Site ID 08-069-0007, located 
~100 km northwest of Denver) exceeded the NAAQS of 75 ppbv as a result of SIs [US 
EPA AQS database, 2017]: March 26th, April 6th, April 27th, May 26-28th, and June 
14th. Several other suburban and rural monitoring stations in Colorado also reported 
exceedences related to SIs on these and other dates in spring 2012 [US EPA AQS database, 


2017]. Less than half of the diurnal variation in the hourly MAMJ O3 at RMNP can be 
explained by the 1st (diurnal) harmonic, therefore other drivers in the O3 variability must 
be considered. Deep SIs, those which impact surface O3 concentrations, were anomalously 
frequent in the western USA in the spring of 2012 compared to the 1990-2012 period [Lin 
et al., 2015], and observed MDA8 O3 was also found to have a maximum in the west- 

ern USA that spring for the period 2004-2012 [Baylon et al., 2016]. This study will ex- 
plore the representation of two of the SIs in remote sensing observations and the God- 
dard Earth Observing System Model, Version 5 (GEOS-5) model and assimilation prod- 
ucts. The first case study in early spring (March 26th local time (LT; +07:00 UTC) here- 
after will be referred to as the SI-1 event and the second case study, which occurred in 
late spring (May 26-28th event LT), will be referred to as the SI-2 event. 


Daily total column O3 (TCO) and relative humidity (RH) from the Atmospheric 
Infrared Sounder (AIRS) on NASA’s Aqua satellite are used to identify the presence of 
SIs over RMNP in observational data and to validate MERRA-2 reanalysis TCO since 
the AIRS O3 data were not assimilated in MERRA-2. AIRS is equipped to measure both 
meteorological variables and chemical profiles [Aumann et al., 2003; Susskind et al., 2006; 
Chahine et al., 2006] and observes the surface twice daily (01:30 and 13:30 LT). The re- 
trievals are performed even when clouds are present which makes the dataset ideal when 
analyzing regions near mid-latitude cyclones. The AIRS team produces several datasets 
of different spatiotemporal resolution. We use the level 3 version 6 (L3 V6) at 1° hor- 
izontal resolution [AIRS Science Team/Joao Texeira, 2013). 


2.2 Model datasets 


NASA’s MERRA-2 reanalysis is an ideal candidate to explore the vertical struc- 
ture of the SIs over RMNP as it is a publicly-available, high-resolution reanalysis dataset 
(0.5° x 0.625° latitude-by-longitude grid, nominally ~50 km in the latitudinal direction, 
72 model layers up to 0.01 hPa [Bosilovich et al., 2016; Gelaro et al., 2017]) which as- 
similates both O3 and meteorological observations [Bosilovich et al., 2015; McCarty et al., 
2016; Gelaro et al., 2017]. The MERRA-2 reanalysis covers the period from January 1, 
1980 to within a couple weeks of real time and is the product of the GEOS-5 data as- 
similation system (DAS) [Bosilovich et al., 2015; Gelaro et al., 2017]. The GEOS-5 model 
includes monthly-averaged ozone production and loss rates linearly interpolated to daily 
values for both the stratosphere and the troposphere [Bosilovich et al., 2016]. After 2004, 
MERRA-2 assimilates satellite retrievals of TCO from the Ozone Monitoring Instrument 
(OMI; Levelt et al. [2006]) and stratospheric O3 profiles from the Microwave Limb Sounder 
(MLS; Waters et al. [2006]) |Bosilovich et al., 2015; McCarty et al., 2016; Gelaro et al., 
2017]. MERRA-2 Og in the lower stratosphere is well-represented and has been shown 
to agree with ozonesondes [Wargan et al., 2015, 2017]; therefore, where there is direct 
influence of stratospheric O3 into the troposphere, such as an SI, we can expect realis- 
tic intrusions although possibly biased since the background ozone in the troposphere 
is simulated by simple chemistry parametrization [Ott et al., 2016]. 


The meteorological and chemical variables — winds (u, v), vertical velocity (w), equiv- 
alent potential temperature (6-; calculated from temperature and specific humidity), Er- 
tel’s PV (EPV), RH and O3 mixing ratios — were extracted on pressure levels up to 150 
hPa [GMAO, 2015a]. In addition, MERRA-2 sea-level pressure (SLP) [GMAO, 2015b] 
and TCO [GMAO, 2015c] are used in the comparison to the AIRS retrievals. 


Two additional model variables are used: GEOS-5 simulated CO using emissions 
described in Ott et al. [2010] and an idealized stratospheric “influence” tracer (STFR) 
from Ott et al. [2016]. The STFR is set to 1 in the stratosphere and to 0 at the surface. 
For the STFR simulation, the tropopause was the higher height of the thermal tropopause 
or the dynamical tropopause. In the GEOS model, the dynamical tropopause is defined 


as the 3 PVU isosurface, where 1 PV unit (PVU) = 10~° K m? kg™! s7!, which is higher 
than the conventional 2 PVU isosurface [Holton et al., 1995]. 


3 Results 


The 8-hour running average Og and the hourly average O3 observed at RMNP and 
the corresponding 3-hourly surface Os from MERRA-2 [GMAO, 2015d] are presented 
here for spring of 2012 (Fig. 1b). During the SI-1 event, the observed hourly O3 at RMNP 
> 75 ppbv for 7 hours (1 hour March 26th and 6 hours March 27th, based on UTC, Fig. 1a) 
with the maximum observed hourly O3 equal to 87 ppbv observed at March 27, 2012 00OUTC 
(Fig. lb; MERRA-2 O3 = 58 ppbv). The second intrusion event, SI-2, influenced ground- 
level Og for several days at RMNP; observed hourly O3 > 75 ppbv for 11 hours during 
this 3-day period (Fig. la) with maximum observed hourly O3 of 91 ppbv at 09 UTC 
on May 27, 2012 (Fig. 1b; MERRA-2 O3 = 65 ppbv). Considering the new NAAQS value 
for MDA8 O3, RMNP observed hourly O3 > 70 ppbv for 11 hours during the SI-1 event 
and 30 hours for the SI-2 event (Fig. la). Figure 1 highlights the doubling of possible 
exceedence days if the new MDA8-O3 NAAQS of > 70 ppbv is applied to 2012. As ex- 
pected, the MERRA-2 surface O3 for the grid box closest to RMNP underestimates the 
O3 variability of a point source measurement (r? = 0.34, based on 968 3-hourly timesteps; 
Fig. 1b) in part because of the simple O3 chemistry in the GEOS-5 model; however, there 
are spikes in the reanalysis O3 at or near the times of observed O3 exceedences, portray- 
ing the influence of stratospheric O3 on the grid-box. 


At the time of an intrusion, relatively dry air is expected to descend toward the 
surface behind a cold front [e.g., Bethan et al., 1998; Cooper et al., 2001; Knowland et al., 
2015]. Due to the topography, SLP over the Rocky Mountains can be difficult to inter- 
pret, however both of the SI case studies occurred when there was a low pressure in the 
Northern Plains region (Fig. 2a,c). During the SI-1 event, there were two low pressure 
systems, one in southwest Montana and one in southeast Wyoming (Fig. 2a). The 700 hPa RH 
was low to the west of a surface trough extending from Wyoming approximately due south 
into Mexico. For the late spring SI-2 event, a cyclone tracked northeastward into North 
Dakota with a cold front trailing into western Kansas where it transitions to a station- 
ary front (Fig. 2c). Here, a new low pressure system formed in southeastern Colorado, 
from which a dry line extends southward through Texas. While relatively low RH is ob- 
served by AIRS to the west of the cold front through the Dakotas and Nebraska, there 
is an even stronger gradient in RH across the dry line (Fig. 2c). 


The SI events can be identified by concurrent observations of O3-rich air with the 
low RH. This can be achieved by focusing on regions where the gradients in TCO are 
large [Olsen et al., 2000; Ott et al., 2016]. The spatial distributions in AIRS TCO and 
MERRA-2 TCO at the approximate time of the AIRS observations agree well [Ott et al., 
2016], although the MERRA-2 TCO is generally biased low compared to observations 
in Fig. 2b,d. This aligns with the findings of Wargan et al. [2017] that MERRA-2 TCO 
in the mid-latitudes was biased low compared to independent TCO measurements from 
the TOMS (Total Ozone Mapping Spectrometer; Herman et al., 1991) instrument. The 
maximum TCO — in both AIRS and MERRA-2 — stretches from the Pacific Northwest 
into the Rocky Mountain states linearly in Fig. 2b and with curvature in Fig. 2d. The 
location of large TCO gradients in Fig. 2b,d correspond to the low RH regions in Fig. 2a,c; 
in particular, large TCO gradients in both AIRS and MERRA-2 and low RH are co-located 
over Colorado (Fig. 2). 


We look for further evidence of the SI-1 and SI-2 tropopause folding events in the 
MERRA-2 reanalysis at the time of maximum O3 at RMNP (Fig. 3). From 300 to 500 hPa 
over the western USA, there are fine-scale filaments of stratospheric air, specifically high 
levels of O3 within the 2 PVU contour, which distinguish the SI events from the back- 
ground (Fig. 3). At the time of maximum hourly O3 observed at RMNP during the SI-1 
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a) Number of hours in an exceedence day (time in UTC) where the RMNP ob- 
served hourly average O3 > 75 ppbv (pink circles) and > 70 ppbv (grey circles) b) 8-hourly 


Figure 1. 


running average O3 (solid black line) and hourly average O3 (dash-dot blue line) from the EPA 
surface observations at RMNP and the 3-hourly MERRA-2 reanalysis surface O3 at the near- 
est grid point to RMNP (40°N, 105.625°W; orange line) for 1 March — 30 June 2012 (time in 
UTC). The exceedence events where the MDA8 O3 > the EPA standard in 2012 (75 ppbv; dot- 
ted horizontal red line) are indicated by the vertical pink shading and the events that would 

be considered exceedences under the new EPA standard (70 ppbv; solid horizontal red line) are 
indicated by the vertical grey shading. The times of the SI-1 and SI-2 events, corresponding to 
Figs. 2-4, are indicated by the black diamonds. 


and SI-2 events, the SI-1 is linear — stretching from Vancouver Island, Canada to the Wyoming- 
Colorado border — as opposed to the curved SI-2 from Washington down to Arizona and 
back to Montana (400 hPa, Fig. 3). Both SI events are a result of a cut-off low (COL) 
near the west coast of the USA in the days prior to the O3 exceedences (not shown). Prior 
to the exceedence at RMNP as a result of SI-1, the tropopause fold rotates around the 
COL and at 500 hPa has a hooked shape off the coast of California (not shown) before 
becoming deformed and elongated — impacting RMNP ~ as the center of vorticity moves 
east and the western portion is being pulled west as a consequence of an Aleutian low 
(Fig. 3a-d). The hook shape of the SI-2 led to the longer period of high O3 (> 75 ppbv) 

at RMNP compared to the duration of high O3 observations associated with SI-1; as the 
SI-2 fold continued to rotate over the western USA at the end of May, there was con- 
tinued draw down of stratospheric air toward the surface over the area, unlike the SI-1 
event which was steered to the northeast as it decayed. It is worth noting that the tro- 
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Figure 2. (a,c) AIRS 700 hPa RH (color; %) with MERRA-2 SLP (4 hPa intervals) 

and (b,d) TCO from AIRS (color; Dobson unit (DU)) and MERRA-2 (dashed; 10 DU 
intervals) for (a,b) SI-1 and (c,d) SI-2. The approximate location of low pressure cen- 

ters (“L”) and frontal boundaries — cold front (line with filled triangles (a,c)), surface 

trough (dashed line (a)), stationary front (filled triangles and half circles on opposite sides 

of line (c)), dry line (line with open half circles (c)) — are presented from the 183UTC 

surface analysis (close to 13:30 LT pass) on (a) March 26, 2012 and (c) May 27, 2012 
(www.wpc.ncep.noass.gov/archives/web_pages/sfc/sfc_archivemaps, Accessed 8 November 2016). 
Note, not all fronts from the analysis archives have been depicted. The location of RMNP (pink 


open circle) and the Colorado state border (think black line) are emphasized. 
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pospheric background levels of MERRA-2 Os are qualitatively consistent with a seasonal 
increase in photochemical production from March (Fig. 3a-d) to May (Fig. 3e-h). 
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Figure 3. Os distribution (color; 5 ppbv increments up to 100 ppbv and increment size in- 
creases above 100 ppbv), geopotential height (thin black contours; 5 dam intervals) and dy- 

namical tropopause (2 PVU isosurface; thick black contour) on 300 (a,e), 400 (b,f), 500 (c,g) 
and 700 hPa surfaces (d,h) corresponding to the time of maximum O3 observations at RMNP 
during the SI-1 event (March 27, 2012 OOUTC; a-d) and the SI-2 event (May 27, 2012 09UTC; 
e-h). Light and dark pink color intervals highlight the previous and current EPA O3 standard, 


respectively. The white dashed lines correspond to transects in Fig. 4. 


The intrusion of air from the stratosphere into the troposphere is captured in ver- 
tical transects for the SI events in the MERRA-2 reanalysis dataset and supported by 
the additional GEOS-5 CO and a fraction of stratospheric air tracer (STFR) data (Fig. 4). 
In the N-S transect through SI-1, the dynamical tropopause and high O3 (> 85 ppbv, 
Fig. dab; > 40% STFR, Fig. 4c) reached altitudes as low as ~600 hPa, and elevated O3 
(> 55 ppbv) reached the surface. Relatively dry air is found within the tropopause fold 
and in the troposphere to the south of the fold (RH < 30 %, Fig. 4b). The low CO (< 110 ppbv, 
Fig. 4d) reached 500 hPa within the tropopause fold; however, since the gradient in CO 
at the base of the fold is less than the gradient of Os, the influence of stratospheric CO 
was lost due to mixing with tropospheric air characterized by higher CO mixing ratios. 
A strong jet at 350 hPa (u > 50 ms~+, Fig. 4e) connects down to the surface. There 
is also a clear frontal boundary on the north-side of the fold extending from an upper- 
level front (indicated by tight isotherms, Fig. 4a-e) down to the surface with strong de- 
scent (w > 80 hPa h~!, white contours, Fig. 4a) and strong ascent (w < -80 hPa h~?, 
black contours, Fig. 4a). 
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Figure 4. Vertical transects of (a-e) the SI-1 event and (f-o) the SI-2 event taken at 

the times of maximum O3 observation at (a-j) 105.625°W from 20°-55°N over RMNP 

(black dot, 40°N) and (k-o) 40°N from 130°-90°W over RMNP (black dot, 105.625°W). O3 
(a,b,f,g,k,l; ppbv), Stratospheric fraction (c,h,m; %), CO (d,i,n; ppbv), Zonal winds (e,j; 
ms‘) and Meridional winds (0; m s~*) are all shown in color with 6. (dashed contour lines, 5 K 
intervals) and the isosurface of 2 PVU (thick black contour). In addition, w (a,f,k; solid contour 
lines, 10 hPa h! intervals, with white contours for descent and black contours for ascent) and 


RH (b,g,l; hatching < 30%) are drawn. Orography indicated by grey region. 


Both N-S and West-East (W-E) transects are shown for the SI-2 case study (Fig. 4f 

o). Although the N-S transect is shown to be just on the eastern edge of the tropopause 
fold (Fig. 3e-h) and the tropopause does not appear to be depressed below 350 hPa over 
RMNP (Fig. 4f-j), there are still strong indicators of a fold in the region. Specifically, 
between 35° to 45°N, there are increased levels of O3 reaching the surface (> 65 ppbv, 
Fig. 4f,g) within an area marked by strong descent (w > 60 hPa h7!, Fig. 4f), low hu- 
midity (RH < 30 % south of 45°N, Fig. 4g), large STFR (> 60 % at 600 hPa, Fig. 4h), 
and low CO (< 110 ppbv, Fig. 4i). The frontal boundary can be identified by the large 
gradients of w and @, to the north of RMNP. This transect highlights the ascent ahead 


of the front (reaching up to ~200 hPa at 50°N, maximum w < -120 hPa h~! at 400 hPa) 
and to a lesser extent the descent behind the front (Fig. 4f). 


Due to the curvature of the SI-2 fold, the W-E transect intersects both sides of the 
hook as seen at 400 hPa in Fig. 3f. The W-E transect captures the tropopause fold at 
105° W over RMNP reaching 550 hPa as well as the western portion of the fold reach- 
ing below 400 hPa at 120°W (Fig. 4k-o). Figure 40 shows the strong jet on both sides 
of the curved tropopause fold (v > 40 m s~! at 105°W and v < -30 ms~! at 120°W). 
The front above RMNP is also seen in this transect by the large gradients in both w (Fig. 4k) 
and 6. (Fig. 4k-o). Specifically, isentropic descent above RMNP brings dry, Oz-rich air 
from the stratosphere towards the surface (O3 > 80 ppbv and STFR > 70 % at 600 hPa, 
O3 > 65 ppbv at surface, Fig. 4k-m). It is interesting to note that low CO (< 110 ppbv, 
Fig. 4n) is simulated to reach the surface at RMNP, despite the large CO values to the 
east emitted by a nearby fire. Biomass burning emissions used in the GEOS-5 simula- 
tion of CO follow the Quick Fire Emission Dataset (QFED) version 2.4r6 which is based 
on MODIS satellite fire radiative power (FRP) [Darmenov and da Silva, 2015]. 


4 Conclusions 


Stratospheric intrusions have been the interest of decades of research, especially 
for the potential influence on ground-level O3 concentrations. However, until recently, 
the fine-scale nature of the O3 filaments have been misrepresented in models and reanal- 
yses, as the features of an SI are best identified in horizontal resolutions of 50 km or smaller 
[Biker et al., 2005; Lin et al., 2012; Ott et al., 2016]. For this reason, and likely because 
reanalysis Og corresponds better with independent observations in the stratosphere than 
in the troposphere [Dragani, 2011; Wargan et al., 2015, 2017], there are very few stud- 
ies of stratosphere-to-troposphere transport which use reanalysis O3 [Skerlak et al., 2014; 
Zanis et al., 2014; Knowland et al., 2015, 2017; Ott et al., 2016; Ryoo et al., 2017]. NASA’s 
MERRA-2 reanalysis is such a high-resolution dataset, which benefits from assimilated 
Os to present O3 on the same spatiotemporal resolution as the meteorology. Here, two 
case study examples of SI events which were known to impact surface O3 air quality are 
examined. The SI events are diagnosed by the folding of the tropopause under the jet 
stream and subsequent isentropic descent of dry, O3-rich/CO-poor stratospheric air to- 
wards the surface using the MERRA-2 reanalysis in combination with surface O3 and 
satellite observations and GEOS-5 simulated CO and a stratospheric tracer. We show 
that MERRA-2, a publicly-available dataset, can be used in scientific studies to iden- 
tify SIs by both atmospheric dynamics and composition. This is a proof of concept study 
opening the door to detailed multi-year analyses of stratospheric intrusions over the USA 
and worldwide. Though the MERRA-2 reanalysis tends to underestimate the magnitude 
of surface O3 during the SIs [see also Ott et al., 2016], the combination of meteorolog- 
ical variables and Og for a relatively long period of time to within a few weeks of present 
time may provide a valuable and unique tool for air quality managers [Kaldunski et al., 
2017] and scientific studies of stratospheric intrusions. 


It is important to be able to identify the differences in anthropogenic and natural 
sources of Os, especially on exceedence days. Since the GEOS-5 model used to produce 
MERRA-2 does not simulate full O3 chemistry in the troposphere, we are unable to de- 
termine the influence of stratospheric O3 on surface concentrations separate from photochemically- 
produced Os, especially in late spring/early summer. The impact of photochemically- 
produced O3 on total Og later in the spring will be explored in more detail using the GEOS-5 
chemistry climate model in a future publication. Yet this study presents strong evidence 
that the MERRA-2 reanalysis can be used in the identification of SIs. 
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